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Abstract 



We investigate the accuracy of several schemes to calculate ground-state 
correlation energies using the generator coordinate technique. Our test-bed 
for the study is the sd interacting boson model, equivalent to a 6-level Lipkin- 
type model. We find that the simplified projection of a triaxial generator 
coordinate state using the S3 subgroup of the rotation group is not very 
accurate in the parameter space of the Hamiltonian of interest. On the other 
hand, a full rotational projection of an axial generator coordinate state gives 
remarkable accuracy. We also discuss the validity of the simplified treatment 
using the extended Gaussian overlap approximation (top-GOA), and show 
that it works reasonably well when the number of boson is four or larger. 
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I. INTRODUCTION 



The self-consistent mean-field theories with a phenomenological nucleon-nucleon interac- 
tion have enjoyed a success in describing ground state properties of a wide range of atomic 
nuclei with only a few adjustable parameters (see Ref. [1] for a recent review). They are now 
at a stage where the ground state correlations beyond the mean-field approximation have to 
be taken into account seriously. This is partly due to the fact that much more accurate calcu- 
lations have been increasingly required in recent years because of the experimental progress 
in the production of nuclei far from the stability line, where the ground state correlation 
beyond the mean-field approximation may play an important role. The major part of the 
correlations produces effects which have smooth trends with proton and neutron number. 
These are already incorporated into the energy functionals of effective mean-field models as, 
e.g., Skyrme-Hartree-Fock or the relativistic mean- field model. However, the correlations 
associated with low-energy modes show strong variations with shell structure, and cannot be 
contained in a smooth energy-density functional. This concerns the low-energy quadrupole 
vibrations and all zero-energy modes associated with symmetry restoration. In fact, the 
correlation effects appear most dramatically for these symmetry modes as there are: the 
center of mass localization, the rotational symmetry, and the particle number conservation. 
Those correlation effects must be taken into account explicietly in order to develop a global 
theory which can be extrapolated to the drip-line regions. 

There are many ways that correlation energies can be calculated. In Ref. [2], we in- 
vestigated a method which uses the random phase approximation (RPA). We found that 
the RPA provides a useful correlation around a spherical as well as for a well deformed 
configurations, but it fails badly around the phase transition point between spherical and 
deformed. Because of this defect, the RPA approach does not seem the best method for a 
global theory. Recently, we have developed an alternative method, called the top-GOA, to 
calculate the ground state correlation energies based on the generator coordinate method 
[3]. This is a generalization of the Gaussian Overlap Approximation (GOA) by taking into 
account properly the topology of the generator coordinate [4]. This method can be easily 
applied to the variation after projection (VAP) scheme, where the energy is minimized after 
the mean-field wave function is projected on to the eigenstates of the symmetry [5]. We 
have tested this method on the three-level Lipkin model, which consists of one vibrational 
degree of freedom and one rotational [3], and have confirmed that the method provides an 
efficient computational means to calculate ground state correlation energies for the full range 
of coupling strengths. 

In this paper, we continue our study on the correlation energies using a model which 
contains the full degrees of freedom of quadrupole motion. To this end, we use a sd inter- 
acting boson model (IBM) [6,7], which may be viewed as a 6-level extension of the Lipkin 
model [8]. The IBM is particularly tailored for the description of the low-lying collective 
modes, thus provides a good testing ground for the present studies of correlations. In real- 
istic systems, treating all the five quadrupole degrees of freedom is a difficult task in many 
aspects. Even if one restricts oneself to the rotational degrees of freedom, one in general has 
to deal with integrals over the three Euler angles, 0, 9, and x- The full triaxial projection 
is still too costly, since a number of rotated wave functions may be required in order to get 
a converged result. Also, the top-GOA scheme for triaxial nuclei is not as simple as in the 
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three-level Lipkin model, because one has to take into account properly the coupling among 
the three Euler angles. How can one overcome these difficulties? We shall study here two 
approximate projection methods. One is the approximate angular momentum projection 
proposed by Bonche et al. [9], which uses the S3 subgroup of the rotation group. With 
this approximation, one needs only five rotated wave functions. The other scheme which we 
consider is the axially symmetric approximation, where the energy is minimized with respect 
to the deformation (3 only, setting the triaxiality 7 equal to zero. With this approximation, 
the integrations for the and \ angles become unnecessary, reducing the projection to a 
one-dimensional integral over 9. The axially symmetric approximation has been widely used 
in the mean-field calculations [10,11], where the approximation seems reasonable given that 
the most nuclei do not have a static triaxial ground state. However, it is not obvious whether 
the approximation remains valid when the fluctuations around the mean-field configuration 
are included, especially when the deformation is small. 

The paper is organized as follows. In Sec. II, we set up the model Hamiltonian and 
discuss several approaches. These include the mean-field approximation, the full triaxial 
angular momentum projection and its approximation, the axially symmetric approximation, 
and the top-GOA for the axial projection. In Sec. Ill, we compare these schemes with the 
exact solutions of the Hamltonian obtained from the matrix diagonalization. We especially 
focus on the feasibility of each method in realistic systems. We then summarize the paper 
in Sec. IV. 



II. SD BOSON HAMILTONIAN 

Consider a N boson system whose Hamiltonian is given by, 

H = H + V = e£4^ - l -Y.QlQ»- (!) 

The first term expresses the single-particle Hamiltonian H , while the second term is the 
residual quadrupole-quadrupole interaction. The quadrupole operator is defined as 

Q M = A 1 ( a td /l + dt a ) + A 2 [dtJ|(^), (2) 

where rf M = (— ) M cL M . When Ai = A 2 = and e > 0, the ground state is the s-boson 
condensed state, whose wave function is given by (s^) N /y/~N\ | ). For a finite value of Ai and 
A2, the Hamiltonian may be diagonalized using the number basis given by, 

\{n}) = \n s n d _ 2 n d _ x n do n dl n d2 ), (3) 

taking only the configurations satisfying 

n s + n d _ 2 + + n do + n dl + n d2 = N, (4) 
-2n d _ 2 - n d _ x + n dl + 2n d2 = 0. (5) 

The first condition, (4), constrains the boson number, while the second equation, (5), is the 
condition that the z component of the angular momentum is zero. With these constraints, 
the basis has a dimenion of 5 for = 2, 18 for N = 4, and 203 for N = 10. We are going 
to compare the exact solutions obtained in this way with results of the collective treatment 
based on the mean-field approximation plus angular momentum projection. 
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A. Mean-field approximation 



We first solve the Hamiltonian in the mean-field approximation. To this end, we consider 
an intrinsic deformed mean-field state given by [7] 

\fr) = : j=(bt) If \), (6) 
where the deformed boson operator is defined as 

The parameter /3 accoutns for the global deformation and 7 for triaxiality. The deformation 
energy surface then reads [7] 

£mf(/3,7) = </W/37>, (8) 
N(3 2 



l + P 2 

+(W-l)|4 i 3 2 - v /|^/J 3 coB3 7 + ?|^| }. (9) 

One finds that the energy minimum appears on the prolate side {j3 > 0,7 = 0) when 
A 2 /Ai < 0, while it is on the oblate side (f3 > 0,7 = 7r/3) for A 2 /Ai > 0. When A 2 is zero, 
the energy surface is independent of 7, corresponding to the 7-unstable case. 



B. Triaxial Angular Momentum Projection 

When (3 is non-zero, the intrinsic wave function (6) is not an eigenstate of the total 
angular momentum J. One can project this state onto the J = state as [5] 

dnR(n)\p 1 )= d<P d X smed9R(<f),9,x)\W, (10) 
Jo Jo Jo 

where R(Q) is the rotation operator. The corresponding energy is given by 

/dSi(/3 7 |B(Si)|/J7> 

Notice that the rotated wave function can be expressed in terms of the rotated boson operator 

as 

\W)=R(Q)\Pi) = ^=(b\ t ) N \), (12) 
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with 



b R = R(npR~\n) 



S t + /5cos 7 E^o^)4 

m 

+^sin 7 X:(^ 2 (fi) + ^_ 2 (Q)) 4 



(13) 



where _D^ m ,(f2) is the Wigner's D function. The overlaps in the projected energy (11) can 
be expressed in terms of commutators such as 

t 6 ' b ^ = TTW ( x + ^ cos2 7 d ° oW + ^ sin2 7 cos(20 + 2x) + ^ W cos(20 " 2x)] 



The results are 



+v / 2/9 2 sin7cos7 4 (^) [cos(2x) + cos(20)]) . 



I(n) = (p 1 \R({l)\f3 1 ) = [b,tf R } N 
V({1) _ {fii\VR{Sl)\fa) 

m) 



eN\l- 



[b,b 



(14) 

(15) 
(16) 



_iV(iV-l) , , 



mm 



Here, we have used the relation 



[i, B N ] = NB"- 1 [A, B] + -N(N - l)B N ' 2 [[A, B],B] + 



(17) 



(18) 



for arbitrary operators A and B. We give an explicit expression for the quadrupole commu- 
tators, [Q m , b R ] and [Q m ,tf], in the Appendix. 

In practice, one can evaluate the integrals in Eq. (11) as follows. First notice that the 
integration intervals for the x an d 4> angles can be reduced from (0, 2tt) to (0, ir) since the 
K quantum number is even for the intrinsic state (6) [12]. Next, because of the reflection 
symmetry of the intrinsic wave function (6) with respect to the z plane, the integration 
range for the 9 angle can be reduced to (0, tt/2). One can then apply the Gauss-Legendre 
quadrature formula to the 9 integral, and the Gauss- Chebyschev formula to the x an d 
integrals [12,13]. One may also try the simpler Simpson formula. We will check the 
convergence of these formulas in the next section. 



C. Approximate Triaxial Projection with Octahedral Group 

Bonchc et al. have considered an approximation to the triaxial angular momentum 
projection (10) based on the octahedral rotation group, that is a group formed from per- 
mutations of the principal axes of inertia [9]. With this representation, the projected wave 
function (10) is approximated as 
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24 

|/3 7 ,J = 0)«^^|/3 7 >, 
i=i 



(19) 



where Si are the 24 elements of the octahedral group. In our case with states even under 
parity, the octahedral group is reduced to S3, the group of permutations of three objects 
(the x,y,z axes). The 6 rotations to be treated are [9], 



Si 


= #(0,0,0) = 1, 


S2 


= £(7r,7r/2,0), 


S3 


= £(-71-/2,-71/2,0), 


S4 


= i2(7r/2,-7r/2,7r/2), 


s 5 


= i?(0,7T,7T/2), 


s 6 


= R(0,ir, -tt/2). 



D. Axial Projection 

When the triaxiality 7 is zero, the and x integrals in Eq. (10) become trivial. The 
triple integral is then reduced to a much simpler single integral with respect to the angle 9. 
This simplifies the projected energy (11) to 

proj ^ J " J 1 _ 1 d(cos9)I(9) ' 1 } 

where the overlaps in this axial approximation read 

^)=( r ^ 27 v( 1 + y( 2 -3-n^)) 7V , (22) 
Ho(9) _ /3 2 (l-fsin 2 fl) 

7(0) l + /3 2 (l-fsin 2 £)' 1 j 

V(0) A 1 



1(9) 2 [1 + /3 2 (1- 



I sin 2 



x 



{ (l + P 2 (l - \ sin 2 9)) (5A 2 + (A 2 + A 2 )/? 2 (l - \ sin 2 (9) 



+(A - l){3 2 ( A 2 (l + 3 cos 2 9) + -!=AiA 2/ 3(l - 3 cos 2 9) 



+^ 2 (4-9sin 2 flcos 2 ^ j. (24) 

The axiallay projected energy (21) depends, of course, only on the global deformation (3. 
VAP means then to minimize the projected energy with respect to the deformation parameter 
(3. 
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E. Top-GOA for Axial Projection 



A further simplification may be achieved using a second-order approach, the extended 
Gaussian Overlap Approximation (top-GOA). In this scheme, the overlaps are expanded 
up to second order derivatives with respect to the generator coordinate while retaining its 
topology. For the axial projection considered in the previous subsection, the procedure is 
very similar as in Ref. [3] for the three-level Lipkin model. From Eqs. (22 - 24), it is clear 
that a natural choice for the expansion variable is sin6>. Expanding the overlaps with respect 
to sin 9, one obtains 

/(e^exp^l^sin 2 ^, (25) 

Mi±m„ EMF(p)+ !m^ e , (26) 

where E MF (P) is the mean-field energy given by Eq. (9) (with 7 = 0), and H 2 (f3) is defined 
as 



d 2 H (6) + V(9) 



d6 2 1(6) 



(27) 



9=0 



Note that we have exponentiated the normalization overlap 1(9) following the idea of the 
Gaussian overlap approximation [14]. 



III. NUMERICAL RESULTS 

A. Comparison of projection schemes 

The exact ground state for the model Hamiltonian (1) and the various integrals needed 
for the projection schemes are solved numerically by standard methods. Figure 1 compares 
the exact solution of the Hamiltonian with the several approximations to the triaxial angular 
momentum projection for N=4 and e = 1. The interaction strength A2 is set to be A2/A1 = 
-V7/A for each Ai, that is a half the SU(3) value, (A 2 /Ai) SU (3) = -V7/2 [15,16]. The 
top panel of the figure shows the ground state correlation energy, i.e., a difference between 
the ground state and the mean-field energies, as a function of the interaction strength Ai. 
The mean-field energy is obtained by minimizing the energy surface (9). The optimum 
deformation parameter (3 thus obtained is shown by the thin solid line in the middle panel. 
One sees the phase transition between the spherical and the deformed configurations at 
Ai = 0.47. The results of full triaxial angular momentum projection, obtained by minimizing 
the projected energy surface (11), are shown by the solid circles in the top panel. These 
results reproduce well the exact results, indicating that the vibrational contribution is not 
large in this model. The optimum deformations (3 and 7 are shown by the thick solid line in 
the middle and the bottom panels, respectively. In contrast to the mean-field approximation, 
the optimum deformation (3 is finite for all the values of Ai, showing no phase transition 
[3]. This is a well-known feature of the variation after projection (VAP) scheme [5]. The 
dotted line in the figure denotes the results of the approximate triaxial angular momentum 
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projection by the S3 subgroup of the rotation group. This method does not seem to provide 
enough correlation energy, and the agreement with the exact results is poor for all the region 



What is the role played by the triaxiality 7 in these calculations? In order to study this, 
we show the results of full axial projection by the dashed line in the figure. These are obtained 
by minimizing the energy function (21), which is equivalent to minimizing (11) while keeping 
7 =0. We find that this approximation reproduces the exact solution remarkably well. The 
result might appear surprising, since the axially symmetric approximation is not expected to 
work near spherical, where all the five quadrupole degrees of freedom should contribute in a 
similar way. However, as we have already discussed, the VAP scheme always leads to a well 
developed deformation even when the mean-field configuration is spherical (see the middle 
panel), and such "dangerous" region can be avoided. Moreover, even though the optimum 
deformation can be small when the interaction strength is very small, this is an irrelevant 
case since the correlation effect is small there. Figure 2 shows the projected energy surface 
Ep TO j(fl, 7), measured with respect to the energy of the pure configuration, s 4 , at A x = 0.5 
and /9=0.741 as a function of triaxiality 7. One sees that the energy gain due to the triaxial 
deformation is indeed small, being consistent with the performance of the axially symmetric 
approximation shown in fig. 1. We summarize the results for Ai = 0.5 in Table I. 

As a further test of the axially symmetric approximation, we repeat the calculations for 
A2/A1 = 0, that is the 7 unstable case. The results are shown in fig. 3, where the meaning 
of each line is the same as in fig. 1. Note that the optimum triaixiality parameter 7 in the 
triaxial angular momentum projection is 30 degree for all the values of Ai, reflecting the 7 
unstable nature of the mean-field approximation. In this case, the performance of the axial 
approximation is not as good as in fig.l (see the dashed line). However, it still provides 
about 80% of correlation energy at Ai =1, and slightly larger at smaller values of Ai, which 
may be acceptable even in realistic systems. 

We notice here that the axially symmetric approximation is sufficient for N = 2 irrespec- 
tive of the value of Ai and A2. From Eqs. (10) and (12), the (normalized) wave function for 
J = state reads 



for any value of 7. The projected wave function is thus independent of 7, and so is the 
projected energy surface. We also note that the axially symmetric approximation becomes 
exact in the limit of N — > 00, as was argued by Kuyucak and Morrison using the 1/N 
expansion technique [17]. For N = 2, the wave function (28) is in fact exact, when (3 is 
minimized. This follows from the observation there are only two J = states in the (sd) A 
configuration space, and their relative amplitudes can be set by a suitable choice of (3, in 
case of attractive interactions. We have checked the trend in between the two limits, N = 2 
and large N. The influence of triaxiality are found strongest around N = 4, where the 
correlation effects are also largest. The effect of triaxiality then decreases slowly as the 
boson number iV increases. 



of A x . 



|/3 7 ,J = 0) 




(28) 
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B. Efficient angular momentum projection 



We next discuss the feasibility of the angular momentum projection. From a computa- 
tional point of view, it is a costly operation to apply the rotation operator to a mean-field 
configuration and take overlaps with it. Thus one wants to minimize the number of points 
in the angular integration mesh. Figure 4 shows the convergence of the angular integrals in 
the projected energy surface (11) with respect to the number of rotated wave functions N TOt , 
for the same parameter set as in fig. 2. Notice that the relations [H, Pj] = and (Pj) 2 = Pj 
are used in deriving Eq. (11), where Pj is the projection operator. For a finite value of 
N TOt , these relations may be violated, and consequently, the numerical formula does not give 
an upper bound of the energy. The open circles are the results of the Simpson method, 
while the closed circles are obtained with the Gaussian quadrature formulas (see Sec. II-C). 
These are for fixed values of deformation parameters f3 and 7, as indicated in the inset of 
the figure. The upper panel is for the axial projection, while the lower panel for the triaxial 
projection. Note that the former is plotted as a function of N Iot , while the latter involves 
the three integrals and is plotted as a function of (N rot Y^ 3 . For the Simpson method, we 
exclude the (0, 6, x) = (0, 0, 0) point in counting the number of state iV ro t in the horizontal 
axis. This state corresponds to the unrotated state from which the rotated wave functions 
are constructed, regardless of which quadrature formula one uses. The figure also shows the 
result of top-GOA and the approximate triaxial projection with the S3 group as a compar- 
ison, which correspond to iV rot = 1 and 5, respectively. From the figure, one observes that 
the convergence for the axial projection is quick if one uses the Gauss-Legendre quadrature 
formula. The energy is almost converged at A" rot = 3. The Simpson method, on the other 
hand, requires more terms to achieve the convergence. For the triaxial projection, a similar 
convergence is seen for each of the three integrals. However, the required number of rotated 
wave functions is as large as 27 in total, making the triaxial angular momentum projection 
with the VAP minimization impractical. The situation is even worse for a larger value of N. 
To demonstrate this, figure 5 shows the results for N = 10. The convergence is somewhat 
slower in this system compared with the N — 4 case. Note that the N rot points-Gauss- 
Legendre formula is exact when the maximum spin in the intrinsic state is J max = 2iV rot — 2 
[12,18]. In the present sd model, the maximum spin J max is given by 2N, and therefore the 
more points are needed in order to get a convergence for the larger value of N. 

Lastly, we discuss the applicability of the top-GOA approach to axial projection (see 
Sec. II-E). This approach requires only one slightly rotated wave function in order to 
evaluate the second derivatives. Figure 6 shows the correlation energy for N = 4 obtained 
with the top-GOA approximation (the dotted line), and with the full axial projection (the 
dashed line). The figure also contains the exact solutions as a comparison. The upper panel 
is for A2/A1 = — a/7/4, while the lower panel is for A2/A1 = 0. We see that the top-GOA 
approximation reproduces the full projection reasonably well. The performance is somewhat 
better for A 2 /Ai = —y/7/4. As was discussed in Ref. [3], the applicability of the top-GOA 
approaches increases quickly for a larger value of boson number N . Indeed, the upper panel 
of figs. 4 and 5 indicates that the agreement between the top-GOA and the exact projection 
significantly improves when A/=10. 
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IV. SUMMARY 



We have used the sci-interacting boson model to investigate projections in a generator 
coordinate approach to calculate the ground state correlation energy associated with the 
quadrupole motion. Our conclusions about the efficiency of various approximations are 
clear. The full angular momentum projection of a triaxial intrinsic state requires a large 
number of rotated wave functions, and it is too costly for realistic calculations. On the 
other hand, we found that the angular momentum projection of an axial intrinsic state 
provides a useful ground state correlation energy. The axially symmetric approximation 
is exact for iV = 2 and iV = oo. The number of rotated wave functions needed there is 
order of 4 if one uses the Gauss- Legendre quadrature formula to compute the angle integral. 
The approximate triaxial projection using the S3 group requires 5 rotated wave functions 
and still performs rather poorly. We thus conclude that the axial projection provides the 
most promising method to compute systematically the ground state correlation energy for 
deformation. 

In applying any projection or generator coordinate expansion, however, one has to bear in 
mind that up to now the energy density functional is defined for a single Slater determinant 
state. It is not designed for a multi-determinantal wave function such as the projected state, 
and there are ambiguities in calculating the density-dependent interaction energy using the 
energy functional. Although several recipes have been proposed, they are all subject to a 
conceptional problem. This difficulty can be avoided in either the following ways. One is 
to use the top-GOA approximation, which can be formulated in terms of the expectation 
values in the mean- field wave function [3]. We have studied the applicability of the top-GOA 
with the present model, and have shown that it already gives a reasonable result for N = 4 
and the performance improves for larger values of N. Alternatively, one may also specify 
the density-dependence in more detail to remove ambiguities. Along these lines, a new form 
of the Skyrme interaction was recently proposed by Duguet and Bonche [19]. In either 
way, the axially symmetric approximation leads to a substantial simplification to perform 
the angular momentum projection with only a few Slater determinants, providing a useful 
means to construct a microscopic global theory for the nuclear binding energy systematics. 
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APPENDIX A: QUADRUPOLE COMMUTATORS 



In this Appendix, we give an explicit expression for the quadrupole commutators, [Q m , b R ] 
and [Q m ,tf], in Eq. (17). For this purpose, it is convenient to use a compact notation for 
the boson operator, b[ m , where b 00 = s and b 2m = d m . Using this notation, we express the 
quadrupole operator Q m and the rotated boson operator b R as 



Qr 



and 



E E mi Ulimz i 



b R = J2B lm (n)bl m , 



(Al) 



(A2) 



respectively. Here, the coefficients ; 2m2 and Bi m (Q) are given by, 

y00,2mi — I J A l S,-m, 

(m) \ r 

?2mi, 00 ~~ A l °mi,mi 

«£l,2m a = (-) ma <2mi2 -m 2 |2m)A 2 , 
5 o(fi) 



From Eqs. (Al) and (A2), one finds 

[Qm,b R ]= E fii^^i^Ci' 

The commutator [Q m , &t] can be obtained by setting Q = in Eq. (A8). This yields, 

[b,QlUQmM] = E E E C^ 2 C^3^ 2 (^)^ 3 (o), 

li,mi h,m2 h,ni3 
b,[Qm,b R ]\ = £ E &,l 2 rnA^MB hmi (0), 

[b,QU,b R ] = E E fi^r^irnxMiWO). 



(A3) 
(A4) 
(A5) 

(A6) 
(A7) 



(A8) 



(A9) 
(A10) 
(All) 
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TABLES 



TABLE I. Comparison of the ground state energy E and the optimum deformation parameters 
(3 and 7 obtained with several methods. The parameters of the Hamiltonian are taken to be N = 4, 
e = 1, Ai = 0.5, and X 2 /\i = -V7/4. The energy is measured with respect to that of the pure 
configuration, s 4 . 



Scheme 


E - E{s A ) 


P 


7 (deg.) 


Exact 


-0.8193 






Triaxial Projection 


-0.8189 


0.741 


17.64 


Axial Projection 


-0.8017 


0.723 


0.0 
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T3 




-1 

-1.5 
-2 
-2.5 



1.6 - 
1.4 1 
1.2 - 
1 

0.8 
0.6 
0.4 
0.2 



\l\= -Sqrt(7)/4 



Exact 

Triaxial Projection 
Axial Projection 
Triaxial Proj. (Approx.) 



Mean Field 
Triaxial Projection 

Axial Projection 

Triaxial Proj. (Approx.) 




FIG. 1. The ground state correlation energy obtained by the several methods (the top panel). 
The parameters of the Hamiltonian are taken to be N = 4, e = 1, and A2/A1 = —y/7/4. The solid 
line is the exact solution of the Hamiltonian obtained by the matrix diagonalization. The dots are 
the results of the full triaxial angular momentum projection, while the dashed line is obtained by 
restricting the intrinsic state to the axially symmetric shape in minimizing the projected energy 
surface. The dotted line denotes the results of the approximate triaxial angular momentum projec- 
tion which uses the S3 subgroup of the octahedral group. The middle and the bottom panels show 
the optimum value of the deformation parameters, (3 and 7, for the angular momentum projections. 
The meaning of the thick solid, the dashed, and the dotted lines is the same as in the top panel, 
while the thin solid line is the result of the mean-field approximation. 
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-60 -40 -20 20 40 60 



Y (deg.) 

FIG. 2. The projected energy surface E pro j((3,^), measured with respect to the energy of the 
pure configuration, s 4 , along the 7 direction for (5 = 0.741. The parameters of the Hamiltonian are 
taken to be N = 4, e = 1, Ai = 0.5, and A 2 /Ai = -V7/A. 



15 




-lh 

-1.5 

-2 
-2.5 h 

-3 

-3.5 
1.4 

1.2 

1 

0.8 

0.6 

0.4 

0.2 h 


40 

35 

30 

25 

20 



Exact 

Triaxial Projection 
Axial Projection 
Triaxial Proj. (Approx.) 

1 1 h 



Mean Field 
Triaxial Projection 
Axial Projection 
Triaxial Proj. (Approx.) 




1 1 

Triaxial Projection 
Triaxial Proj. (Approx.) 







0.5 



FIG. 3. Same as fig.l, but for A 2 /Ai = 0. 
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Triaxial Projection N = 4 




= 0.741, 7=0.308 
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•— • Gaussian Quadrature 






▲ S 3 subgroup 




: \ 






,1,1,1,1,1,1,1,1,1, 



2 4 6 8 10 12 14 16 18 20 

(N rol ) 1/3 

FIG. 4. Influence of the generator coordinate truncation on the ground state energy. The upper 
and the lower panels are for the axial and the triaxial projections, respectively. The former plots 
the energy as a function of the number of rotated wave functions iV ro t , while the latter plots as a 
function of (iVrot) 1 / 3 ; for the optimum values of the deformation parameters (3 and 7 indicated in 
the insets. The open and the closed circles are the results of the Simpson method and the Gaussian 
quadrature formula, respectively. The triangles denote the result of the top-GOA approximation 
(in the upper panel) and that of the approximate projection with the S3 subgroup (in the lower 
panel). The parameters of the Hamiltonian are the same as in fig. 2. 
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FIG. 5. Same as fig.4, but for N = 10. 
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FIG. 6. The correlation energy obtained in the axially symmetric approximation as a function 
of Ai for N = 4. The upper and the lower panels are for A2/A1 = — V7/4 and A2 = 0, respectively. 
The dashed line is the result of the full axial projection, while the dotted line is obtained in the 
top-GOA approximation to the axial projection. The exact solution of the Hamiltonian is denoted 
by the solid line. 
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